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Magnetosphere-Ionosphere Coupling Through E'-region 
Turbulence: Energy Budget 

Y. S. Dimant and M. M. Oppenheim 

Abstract. 

During periods of intense geomagnetic activity, strong electric fields and currents pen- 
etrate from the magnetosphere into high-latitude ionosphere where they dissipate energy, 
form electrojets, and excite plasma instabilities in the B-rcgion ionosphere. These insta- 
bilities give rise to plasma turbulence which induces non-linear currents and strong anoma- 
lous electron heating (AEH) as observed by radars. These two effects can increase the 
global ionospheric conductances. This paper analyzes the energy budget in the electro- 
jet, while the companion paper applies this analysis to develop a model of anomalous 
conductivity and frictional heating useful in large-scale simulations and models of the 
geospace environment. Employing first principles, this paper proves for the general case 
an earlier conjecture that the source of energy for plasma turbulence and anomalous heat- 
ing equals the work by external field on the non- linear current. Using a two-fluid model 
of an arbitrarily magnetized plasma and the quasilinear approximation, this paper de- 
scribes the energy conversion process, calculates the partial sources of anomalous heat- 
ing, and reconciles the apparent contradiction between the inherently 2-D non-linear cur- 
rent and the 3-D nature of AEH. 



1. Introduction 

At high latitudes, the large-scale electric field, Bo, from 
the solar wind and Earth's magnetosphere penetrates to the 
ionosphere. Across the B/upperB-region altitudes, roughly 
between 80 and 130 km, the electrons are strongly mag- 
netized, while the ions become at least partially demagne- 
tized due to frequent collisions with the neutral atmosphere. 
This demagnetization slows the drifting ions down, mak- 
ing the Bo x Bo drift of electrons into strong Hall currents 
named high-latitude electrojets. In the global picture of 
magnetosphere-ionosphere (MI) coupling, this is the region 
where field-aligned magnetospheric currents close and dissi- 
pate energy. 

Electrojet currents often drive plasma instabilities that 
generate electrostatic field fluctuations coupled to plasma 
density irregularities. These fluctuations have low frequen- 
cies, usually smaller than the average frequencies of electron 
and ion collisions with neutrals, f e ,i, while the correspond- 
ing wavelengths in the unstable range exceed the ion mean 
free path. Density irregularities, usually in the range from 
tens of centimeters to tens of meters, are routinely detected 
as strong coherent radar echoes [e.g., Cohen and Bowles, 
1967; Balsley and Farley, 1971; Crochet et al, 1979; Kudeki 
et al, 1987]. Rocket flights through the lower ionosphere 
have detected also electrostatic field fluctuations [e.g., Pfaff 
et al, 1987, 1992, 1997; Rose et al, 1992; Fukao et al, 1998]. 
These irregularities and fluctuations are caused by variety 
of instabilities including the Farley-Buneman (FB) [Farley, 
1963; Buneman, 1963], gradient drift [Hoh, 1963; Maeda 
et al, 1963], and thermal instabilities [Dimant and Sudan, 
1995, 1997; Kagan and Kelley, 2000; Dimant and Oppen- 
heim, 2004]. The FB instability is excited when the relative 
velocity between the average electron and ion streams ex- 
ceeds the local ion-acoustic speed. At high latitudes, this 
usually occurs when Bo = Boj exceeds the threshold value 
of about 20 mV/m. These and much stronger fields are not 
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uncommon in the sub-auroral, auroral, and polar cap areas, 
especially during magnetospheric storms and substorms. 

Small-scale fluctuations generated by these instabilities 
can cause enormous anomalous electron heating (AEH). For 
about thirty years, radars have observed strong electron 
temperature elevations from 300-500 K up to more than 
4000 K, correlating with strong convection fields Eo [Schlegel 
and St. -Maurice, 1981; Providakes et al, 1988; Stauning and 
Olesen, 1989; St.-Maurice et al, 1990; Williams et al, 1992; 
Foster and Erickson, 2000; Bahcivan, 2007]. Simple esti- 
mates show that regular ohmic heating by Eo alone cannot 
account in full measure for such huge temperature eleva- 
tions. A strong correlation between AEH and E , as well 
as other physical arguments, shows that average heating 
by FB-generated turbulent electric fields causes AEH [St.- 
Maurice and Laher, 1985; Robinson, 1986, 1992; Providakes 
et al, 1988; St.-Maurice, 1987, 1990; Dimant and Milikh, 
2003]. AEH occurs largely because the turbulent electro- 
static field, SE = — V5$, has a small component 5Eo par- 
allel to the geomagnetic field Bo [St.-Maurice and Laher, 
1985; Providakes et al, 1988; Dimant and Milikh, 2003; Mi- 
likh and Dimant, 2003; Bahcivan et al, 2006]. The impor- 
tance of (5B|| makes the entire process fully 3-D. 

Anomalous electron heating can modify ionospheric con- 
ductances and hence affect the coupling between the mag- 
netosphere and ionosphere. Any electron heating directly 
affects the temperature-dependent electron-neutral collision 
frequency and, hence, the electron part of the Pedersen con- 
ductivity. This part, however, is usually small compared 
to the electron Hall and ion Pedersen conductivities. How- 
ever, AEH causes a gradual elevation of the mean plasma 
density within the anomalously heated regions by reduc- 
ing the local plasma recombination rate [Gurevich, 1978; 
St.-Maurice, 1990; Dimant and Milikh, 2003; Milikh et al, 
2006] . The AEH-induced plasma density elevations increase 
all conductivities in proportion. However, this mechanism 
requires tens of seconds or even minutes because of the 
slow development of the ionization-recombination equilib- 
rium. If Bo changes faster than the characteristic recom- 
bination timescale then its time-averaged effect on density 
will be smoothed and reduced. 

Plasma turbulence, however, can directly modify local 
ionospheric conductivities via a wave-induced non-linear 
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current (NC) associated with plasma density irregularities 
[Rogister and Jamm, 1975; Oppenheim, 1997; Buchert et al., 
2006] . The physical nature of the NC is explained in Fig. 1 
for magnetized electrons and unmagnetized ions. The wave 
field, SE, has different signs in the wave density maxima 
and minima, so that the corresponding SE x Bo-drifts of 
magnetized electrons have opposite directions. As a result, 
more electrons drift along the maxima than in the opposite 
direction along the minima, producing a net non-linear cur- 
rent density, jsc- At higher B-region altitudes, partially 
magnetized ions can also contribute to this current. Rela- 
tive density perturbations in saturated FB turbulence may 
reach at most tens percent, while the rms turbulent field, 
{SE 2 } 1 ^ 2 , is comparable to Eo (the angular brackets here 
and below denote spatial-temporal averaging). As a result, 
the total NC, considered as a plasma response to the exter- 
nal electric field Eo, amounts to only a fraction of the regular 
electrojet Hall current. However, in most of the electrojet 
the NC is directed largely parallel to Eq, so that it may in- 
crease significantly the much smaller Pedersen conductivity. 
This is critically important because the Pedersen conductiv- 
ity allows the MI field-aligned currents to close and dissipate 
energy. The combined effect of the NC and AEH makes 
the ionosphere less resistive. The anomalous conductance 
may, at least partially, account for systematic overestimates 
of the total cross-polar cap potential in global MHD mod- 
els that employ laminar conductivities [e.g., Winglee et al., 
1997; Raeder et al, 1998, 2001; Siscoe et al, 2002; Ober 
et al, 2003; Merkin et al, 2005a, b, 2007; Guild et al, 2008; 
Wang et al, 2008]. 

Furthermore, the NC and AEH are intrinsically related. 
As discussed in this paper, the work by the external DC elec- 
tric field on the total current equals the total field energy 
input to ionosphere (total Joule heating), including that re- 
sponsible for AEH. On the other hand, the total frictional 




Figure 1. Formation of a net non- linear current (NC) at 
a given wave of plasma compression/decompression with 
the wavevector k in the Eo x Bo-direction. The wave 
electrostatic field, <5B||fe, has opposite directions in the 
plasma density maxima (8n > 0) and minima (8n < 0), 
resulting in the oppositely directed 8E x Bo drifts of 
magnetized electrons. More negatively charged particles 
move in the —Bo than those in the opposite direction, 
resulting in formation of the net positive current, jwc, 
parallel to Bo. 



heating of electrons and ions effectively increases wave dis- 
sipation and reduces the turbulence intensity, thus affecting 
the NC. 

The process of anomalous energy deposition from the 
magnetosphere to ionosphere has been recently studied by 
Buchert et al. [2006]. They showed that the apparent aver- 
age turbulent energy deposition per unit volume and time, 
{SE ■ 8j) , equals zero (here Sj is the fluctuation of the total 
current density j; a number of notations in this paper differ 
from those of Buchert et al. [2006]). The authors deduced 
that the total energy per unit volume and time lost by the 
external field Bo to particles equals the total work performed 
by this field on the total average current, j = jo + Jnc (jo is 
the laminar current density). They suggested that the actual 
energy input for the instability development and turbulence- 
induced energy losses equals Bo • jnc- Their calculations 
supported this remarkable conjecture, but provided no com- 
plete proof since their quantitative analysis was based on an 
oversimplified and restrictive model. First and most impor- 
tantly, they studied only 2-D turbulence in a perpendicular 
to Bo plane. The authors did not state the 2-D restric- 
tion explicitly, but it becomes evident from their Eqs. (17) 
and (18) which effectively exclude ku. Due to this, their 
treatment leaves out the key 3-D energy conversion process 
primarily responsible for AEH. Secondly, they employed the 
simplest two-fluid model with fully demagnetized ions. This 
widely used approximation is good for most of the electrojet 
but fails at its top part, right where the Pedersen conductiv- 
ity reaches its maximum. Thirdly, their calculations of NC 
and other average terms were fully based on a quasilinear, 
narrowband approximation for FB waves. In this approx- 
imation, turbulence consists of relatively small-amplitude 
waves which are described by linear relations with narrow- 
band wave frequencies. In many cases, this is a reason- 
able approach, although its applicability is less justified for 
the driving electric field well above the instability threshold, 
when the effect of AEH is especially strong. 

These restrictions and some inconclusiveness associated 
with them raise a number of important questions. First 
and foremost, is the conjecture about Bo • Jnc as the energy 
input for plasma turbulence really true? If so, then this 
fundamental fact should follow directly from the first prin- 
ciples and be universally applicable. In particular, it has 
to be valid for arbitrarily magnetized particles and strongly 
non-linear processes. Then, is it possible to deduce this 
fact from a general viewpoint with no or minimal approx- 
imations? Second, presuming that Bo ■ Jnc is the correct 
energy input, then how is it compatible with the 3-D effect 
of AEH largely associated with small SEnl This question 
arises because j'nc, as well as Bo, lies almost precisely in 
the plane perpendicular to Bo with no average contribution 
from k\\ . Recently, we performed a number of 2-D and 3-D 
particle-in-cell (PIC) simulations in big periodic boxes with 
unprecedentedly dense meshes and large numbers of PIC 
particles [Oppenheim et al., 2011]. Comparing the 2-D and 
3-D results for the same background parameters, we have ob- 
tained compelling evidence that anomalous electron heating 
related to plasma structuring along Bo does exist. Then, 
how could essentially 3-D heating originate from the 2-D 
energy input? Further, how this energy input is distributed 
between various groups of particles? This is important for 
making accurate estimates of anomalous heating for differ- 
ent particle groups. During strong magnetospheric events, 
what feedback of developed B-region turbulence on global 
MHD behavior of the magnetosphere might be expected? 
How to quantify this effect for including it in global MHD 
codes intended for space weather predictions? Finally, what 
channels provide the corresponding energy flow between the 
magnetosphere and B/B-region ionosphere? 
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In this paper, we address these issues and create a rig- 
orous basis for calculating anomalous conductivities in the 
companion paper [Dimant and Oppenheim, 2011]. We start 
by confirming from first principles that fully saturated tur- 
bulence does yield {5E ■ Sj) = 0, then confirm the Buchert 
et al. [2006] deduction regarding the turbulent energy input 
and establish its universal validity. In order to quantitatively 
develop a 3-D model of AEH and resolve the apparent con- 
tradiction between this interpretation and the 2-D nature 
of the energy input, we perform specific calculations for the 
case of arbitrary particle magnetization, using a quasilinear 
approximation. We calculate the non-linear current, total 
energy input, and partial average frictional heating sources 
for both electrons and ions in terms of a given spectrum 
of density irregularities. We show that the major quantita- 
tive difference between 2-D and 3-D developed turbulence 
lies in the magnitude of density perturbations. These per- 
turbations and the non-linear current proportional to them 
are noticeably larger in 3-D than in 2-D. This difference ex- 
plains the larger energy input in 3-D and is responsible for 
AEH caused by turbulent fields, thus resolving the above- 
mentioned contradiction. 

The paper is organized as follows. In Sect. 2, using only 
first principles with no approximations, we confirm for the 
general case the Buchert et al. [2006] findings regarding the 
energy input. In Sect. 3, we develop a quasilinear approach, 
similar to that of Buchert et al. [2006], but for the general 
3-D case of arbitrarily magnetized particles. This allows us 
to calculates partial non-linear currents, relevant energy in- 
puts, and frictional heating sources. In Sect. 4, we discuss 
global energy flow between the magnetosphere and iono- 
sphere. In the appendix, we check the validity of the conven- 
tional electrostatic approximation for lower-ionosphere wave 
processes. 



2. Energy Conversion: 
Consideration 



First Principle 



In this section, we derive general relations regarding the 
average energy input in quasi-periodic systems and show 
how the Buchert et al. [2006] deductions follow from funda- 
mental electrodynamic and plasma kinetics principles with 
no approximations like electrostatics, quasi-neutrality, fluid- 
model description, etc. 

First, we consider the evolution of the field energy in plas- 
mas by considering the exact electrodynamics, starting with 
Ampere's and Faraday's laws, 



dtE = c V x B 



d t B 



-V x E, 



(la) 
(lb) 



where eo is the permittivity of free space, c is the speed of 
light in vacuum, E and B are the electric field and magnetic 
induction, and j is the total current density. 

Taking scalar products of Eq. (la) with EoE, Eq. (lb) 
with £oc 2 B and adding the results, we obtain the standard 
energy balance equation (aka Poynting's theorem): 



dtU + V-S = -E- j, 



where 



U = e ± [E 1 + c 2 B 2 



E x B 
Mo 



(2) 



(3) 



are the field energy density and the corresponding flux (the 
Poynting vector), respectively; fig = (eoc 2 ) -1 is the perme- 
ability of free space. 

Now we need to find j from the plasma. The dynamics 
of individual particles of type s, such as electrons or ions 
(s = e, i), is accurately described by Boltzmann's kinetic 



equation, which can be written in the 6-D divergence form 

as 



dtfs + V • (v s fs) + dv 3 



-(E + v 3 x B)f s 



= S 3 



(4) 



Here v s is the kinetic velocity of the particles-s, while 
f s (f,t,v s ) is their single-particle velocity distribution func- 
tion normalized to the particle-s density, n s = f f s d 3 v s (the 
integration here and below is performed over the entire 3- 
D velocity space); q s and m s are the particle charge and 
mass; S s is the collisional operator which includes particle- 
s collisions and can also include the ionization sources and 
recombination losses; for simplicity, we disregard effects of 
gravity. Multiplying Eq. (4) by m a v 2 /2, integrating over 
the velocity space, and adding the results for all plasma 
particles, we obtain the energy balance relation for plasma, 



d t ^£ s + V • ^K s = E ■ j + ^L s 



(5) 



Here £ s = f(m s v 2 /2) f s d 3 v s and K s = f(m s v 2 /2)v s f s d 3 v s 
are the average particle energy and energy-flux densities, re- 
spectively; L a = J(m s v 2 /2)S s d 3 v s combines all collisional 
energy gains and losses; j = q s n s V s , is the total current 
density, the same as in Eq. (2); V s = f v a fsd 3 v a /n a is the 
particle-s mean fluid velocity. In the general case, the par- 
ticle energy density £ s combines the mean thermal energy, 
3n s T s /2, where T s = f[m s (v s - V s ) 2 /3]f s d 3 v s /n s is the ef- 
fective particle temperature, with the kinetic energy density 
of the mean particle flow, n s m s V 2 /2. 

Comparison of Eqs. (2) and (5) shows that E-j is the total 
energy input from the fields to particles per unit volume and 
time. The energy deposited in an individual group of par- 
ticles at a given location may be then slightly redistributed 
via Coulomb collisions and transported to other locations. 
Eventually, this energy becomes lost to the abundant neu- 
tral atmosphere via predominantly inelastic plasma-neutral 
collisions. 

General Eqs. (2) and (5) apply to all plasma processes, 
linear or non-linear. To specifically discuss turbulent pro- 
cesses, we need to separate them from the slowly evolving 
and large-scale regular background structures and processes. 
This can be done by using a conventional two-scale approach 
in which macroscopic background structures and processes 
are presumed to have much longer characteristic spatial and 
temporal scales than does the turbulence. The physical con- 
ditions in the i5-region plasma make this procedure more 
applicable to plasma irregularities generated by local insta- 
bility mechanisms, such as the Farley-Buneman (FB) and 
thermal-driven instabilities [Dimant and Sudan, 1997; Ka- 
gan and Kelley, 2000; Dimant and Oppenheim, 2004] than 
for relatively large-scale irregularities generated by plasma 
gradients. 

We assume that plasma turbulence generated by insta- 
bilities consists mostly of waves whose wavelengths \i in 
each direction i are much smaller than the corresponding 
typical scales of spatial variation of the macroscopic back- 
ground parameters, A,. Then, for any given location r, we 
can build an imaginary rectilinear box centered around f 
with the sizes Li that satisfy the conditions \i <§C Li <§C A;. 
The characteristic scales in different directions i can differ 
dramatically. For example, the wavelengths of nearly field- 
aligned irregularities parallel to Bo are generally orders of 
magnitude larger that those perpendicular to Bo, so that 
the corresponding box sizes would also be vastly different. 
Within the imaginary box, we can extend the background 
parameters from the box center r uniformly to the entire box 
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and impose periodic boundary conditions. This approach is 
widely employed in computer simulations. 

We can similarly describe the temporal evolution of the 
process by presuming long-lasting non-linearly saturated 
turbulence in a quasi-stationary background. This implies a 
continuous driving of instability and allows one to introduce 
a quasi-period T, analogous to Li, which is much longer 
than typical timescales of turbulent variations but is much 
shorter than a characteristic timescale of slow evolution of 
the macroscopic background. 

To implement the two-scale procedure, we introduce the 
following convention. We denote the short-scale and fast 
periodic variables for turbulent processes by x = {x\, xi, xz) 
and t, while denoting the corresponding coordinates and 
time of large-scale and slow background variations by r and 
t, respectively. With respect to the periodic coordinates x 
and time t, we define a spatial-temporal average as 



the result in accord with the definition of Eq. (6), we obtain 



(•••>= r r r Tf, / dx i / dx 'i / dx 3 / (• • • ) dt. 

L 1 L 2 L;il J-L 1 /2 J-L 2 /2 J-L3/2 Jt-T/2 

(6) 

In general, such averages can remain functions of the slow 
variables r and r. This description of turbulence does not 
require the ensemble averaging employed by Bucheri et al. 
[2006]. 

Having implemented this procedure, we split the quasi- 
periodic fields and currents into their r,r-dependent average 
parts and the corresponding local periodic perturbations, 

E={E) + SE, B = {B)+SB, j = (J) + Sj, (7) 

Using (SE) = {Sj) — 0, we obtain for the average energy 
deposition term 



{E-j) = {E)-{j) + {5E-5j). 



(8) 



One might naively interpret the first term on the right-hand 
side (RHS) of Eq. (8) as the energy input from the back- 
ground fields to the undisturbed plasma, while the second 
term as the corresponding average contribution from the 
turbulent fields. As mentioned in the Introduction, Buchert 
et al. [2006] found, using a simplified model, that {SE ■ Sj) 
equals zero by the following formal mathematical reason. If 
one expands this term further in terms of the mean particle 
fluid velocities to the lowest-order quadratic non-linearity 
then the expected total turbulent frictional heating term, 
S s q s n a0 {5E ■ SV S ), turns out to be automatically canceled 
by a density perturbation term, Yl a Qs{Sn s 5E) ■ V s o (here 
8n s and SV S are the perturbations of the densities and par- 
ticle fluid velocities). Right below we show that for purely 
periodic and spatially homogeneous turbulence the equality 
{SE ■ Sj) — is a natural and universal constraint, required 
merely by the imposed periodicity. This constraint follows 
directly from exact Maxwell's equations without invoking 
specific plasma models. 

Indeed, Maxwell's equations are linear, allowing the sep- 
aration of the average quantities from wave perturbations. 
The perturbations depend on the small-scale, quasi-periodic 
variables, x, t, and can also have an adiabatically slow r, 
r-dependence due to background inhomogeneities and evo- 
lution. In principle, the inhomogeneous background param- 
eters may evolve in such a way that an instability threshold 
is crossed, resulting in sudden onset or disappearance of in- 
stability in some locations at certain moments of time. Such 
instances, analogous to second-order phase transitions [Lan- 
dau and Lifshitz, 2000], break the validity of our two-scale 
approximation and deserve a special treatment that lies be- 
yond the framework of this paper. Apart from these special 
occasions, we can apply to the perturbations of the fields 
and current the same steps that lead to Eq. (2). Averaging 



e {(SE 2 )+c 2 (5B 2 )) 



{SE x SB) 



Mo 



+ Vy - = -{SE-5j). 

(9) 

Then for saturated turbulence with constant average charac- 
teristics in a strictly periodic box, the left-hand side (LHS) 
of Eq. (9) disappears, yielding 



{SE ■ Sj) = 0, 



(10) 



regardless of the specific fluid or kinetic models employed 
for the plasma description. For strictly periodic processes, 
this result is exact. One can also obtain it directly from the 
discrete Fourier harmonics of the electric field and current. 
We do this right below and also demonstrate that electro- 
static and quasi-neutral approximations do not lift this exact 
electrodynamic constraint. 

We start by introducing spatial and temporal Fourier 
harmonics of wave perturbations for general, linear or non- 
linear, periodic processes. Since we imply saturated tur- 
bulence in a 4-D box (3-D space + time) with periodic 
boundary conditions, it is logical to employ the 4-D dis- 
crete Fourier transformation. For easy reading, we assign to 
the Fourier transforms the notations of the original variables 
but with additional subscripts k, u>. For a scalar or vector 
periodic perturbation 5F(x, t) as a function of coordinates 
Xi and time t, we define such transformations within the 4-D 
box of the corresponding sizes L t and the quasi-period T as 



5F(x,t)= J2 <5^exp[i(fc-f-^)], 

^L 3 f f sF & *) exp ~m kiXi ~ ut 



SF *+> L 



(11a) 



d 3 Xidt. 



(lib) 



Here the k,uj 7^ summation is taken over all discrete val- 
ues of k, ki = 2nni/Li, lo = 2irm/T, with integer m, m, 
while the integration is performed over the entire 4-D box, 
as in Eq. (6). Since we apply the discrete Fourier transforms 
only to the perturbations, we exclude from the summation 
over harmonics the average values corresponding to k, u> = 0. 
In terms of their Fourier harmonics, 8 At and SBz , the 
spatial-temporal average of any product, scalar or vector, of 
two functions A(x,t) and B[x,i), is given by 

{SA(x,t)6B(x,t)) = £ SA i JBl u =Y t SA^JB^. 

(12) 

Our current discrete-transform normalization differs from 
the continuous one in Buchert et al. [2006] and allows to 
avoid the emergence of extraneous factors like VT {V = 
LiL 2 L 3 ) and (2tt) 4 in the explicit physical expressions for 
average quadratically non-linear quantities [e.g., Buchert 
et al, 2006, Eqs. (6)-(8), (19) -(26), etc.]. 

Using Fourier transforms, we now prove Eq. (10) in a 
more direct way. Applying Eq. (11) to Eq. (1), for a given 
Fourier harmonic, we obtain 

-, _ 

ic kxSB K = — 5j K -iuSE K ikxSEj; = iu>8B K 

£0 

(13) 

Making a cross-product of Eq. (13) with k and using 
k ■ SBj: u — 0, we express the field perturbations in terms of 
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ik x Sj^ L 



k,uj 



e cj(fc 2 c 2 — a; 2 ) 



(14) 
(15) 



In a long-lived quasi-periodic non-linearly saturated state, 
all wave frequencies w must be real, so that Eq. (15) yields 
Re(*^ >w • SJIJ = 0, i.e., Eq. (10). 

General Eq. (10) represents a strict constraint that fol- 
lows from the full electrodynamics of quasi-periodic pro- 
cesses, but it is not obvious that it should hold for the elec- 
trostatic and quasi-neutrality approximations. The follow- 
ing demonstrates that these approximations do not lift this 
constraint. 

Indeed, the electron and ion continuity equations com- 
bined yield edt(Sn e — Siii) — VSj, so that ew{Sn i ^ w — 

& n eku) — k- Sjku- Combining this with Poisson's equa- 
tion, £qV • SE = e(Sn e — Sni), we obtain 



ik ■ SE^ 



k ■ 6 k,< 



(16) 



For an electrostatic field, we have 5E^ w — \SE { g u \ (k/k), so 
that Eqs. (12) and (16) yield 

{Sj- SE) = Re(6jr ■ 5EU = -s \5Er J 2 Imw = 0. (17) 



If we add quasi-neutrality, V • Sj = 0, i.e., k ■ 6j% 







for individual harmonics, then the satisfaction of Eq. (10) 
becomes obvious even without assumption of real u, 



{5j-5E)=i J2 (*-tfs,J**E, w = ' 



(18) 



k,u>^0 



Equation (10) does not mean, however, that developed 
plasma turbulence makes no contribution to the average par- 
ticle heating. This would certainly contradict both observa- 
tions and PIC simulations. The paradox can be resolved 
as follows. The average electric field is merely the external 
field, (E) = E , while (j) j , where Jo = J] s q s n a0 V s0 is 
the undisturbed current density determined by the laminar 
plasma response to Eo- In addition to jo, the total average 
current density, (j) , includes a wave- induced direct NC [Ro- 
gister and Jamin, 1975; Oppenheim, 1997; Buchert et al., 
2006], 

T° = 0)-Jo = ^q s {Sn s SV a )- (19) 

s 

The total average loss of field energy to particles per unit vol- 
ume and time is given by {E ■ j). The corresponding energy 
dissipation in the laminar ionosphere with no instabilities is 
given by Eo ■ jo- Then the total energy dissipation exclu- 
sively due to plasma turbulence is Lturb = (E ■ j) — Eq ■ jo- 
Using Eqs. (10) and (19), one can easily establish that 



iturb = P nc = Eo ■ jr-JC • 



(20) 



Thus, it is formally the work by the external electric field on 
the non-linear current, Pnc, rather than {SE • Sj), that pro- 
vides the required turbulent energy deposition for all kinds 
of anomalous heating of plasma particles. Buchert et al. 
[2006] showed this for the restricted case of fluid plasmas 
with fully unmagnetized ions and a quasilinear wave descrip- 
tion, but did not establish it in the general case. We have 



just demonstrated that this fundamental result follows di- 
rectly from the general field electrodynamics and no specific 
plasma models. By their physical meaning and according to 
Eq. (20) each of the two equal quantities L tur b and Pnc can 
be named "turbulent Joule heating." 

Rigorously speaking, exact Eqs. (10) and (20) apply only 
to a homogeneous and stationary background. For mildly 
inhomogeneous and slowly evolving background parameters, 
Eq. (9) can be treated using a regular perturbation tech- 
nique. The zero-order approximation, corresponding to a 
given turbulence level within an isolated box with a uniform 
background and periodic boundary conditions and, hence, 
not affected by the outside inhomogeneity, yields Eq. (10). 
To reach a next-order accuracy, one has to establish the 
zero-order parameter dependence of non-linearly saturated 
turbulence characteristics, using, e.g., a series of computer 
simulations with periodic boundary conditions but various 
background parameters. Given the large-scale spatial de- 
pendence and slow evolution of the background parameters, 
one could calculate then the LHS of Eq. (9). This would 
yield the first-order, potentially non-zero, values for {SE-Sj) 
on the RHS. Under applicability of our two-scale approach, 
however, such possible finite values of {SE ■ Sj) will auto- 
matically be small compared to the nearly balancing each 
other leading terms in the expanded form of {SE ■ Sj) with 
Sj — ~}2 3 q s (V s oSn + n SV s + Sn8V s ), creating only a rela- 
tively small mismatch between L t urb and Pnc. 

Below we calculate energies deposited among individual 
particle groups using an approach similar to that used by 
Buchert et al. [2006], except that we will apply it to a system 
with arbitrarily magnetized electrons and ions. This applies 
to all altitudes across the entire E region down to the upper 
D region. As mentioned above, Buchert et al. [2006] actu- 
ally performed a 2-D treatment. We consider here fully 3-D 
turbulence, which is crucial for anomalous electron heating. 
Note that the total energy input spent on anomalous heat- 
ing of plasma particles is less than Pnc because a fraction of 
the deposited energy via collisions goes directly to colliding 
neutrals and has no chance to heat the plasma. 

3. Partial Energy Deposit: Quasilinear 
Approximation 

In this section, we employ Fourier harmonics in a 3-D 
periodic box in order to identify the effect of the paral- 
lel turbulent electric field. Our results apply to all regions 
that contribute to the total ionospheric conductances, from 
the top electrojet down to the potentially unstable Z)-region 
[Kelley, 2009]. Similar processes can occur in other plasma 
media, like the Solar chromosphere [Liperovsky et al, 2000; 
Fontenla et al, 2008; Gogoberidze et al, 2009], other plan- 
etary ionospheres, and laboratory plasma [D'Angelo et al., 
1974; John and Saxena, 1975; Koepke, 2008], despite the 
dramatic differences in parameters. Lastly, these calcula- 
tions for arbitrarily magnetized plasma serve as an addi- 
tional verification of the general relations obtained in the 
previous section from first principles. 

In this paper, we restrict ionospheric particles to electrons 
and a single species of ions, adopting the quasi-neutrality, 
n e ~ n, = n. Similarly to Buchert et al. [2006], we apply 
here a quasilinear approximation, by which we imply that 
separate Fourier harmonics of predominantly electrostatic 
field and plasma density fluctuations are coupled through 
simple linear relations. Using these relations, we calcu- 
late then quadratically non-linear averages in terms of given 
spatial-temporal turbulence spectra. 

This section is organized as follows. In Sect. 3.1, we ob- 
tain linear relationships to first-order accuracy and explain 
what we mean by the different orders. These first-order re- 
lations provide the linear wave frequencies and relations be- 
tween the electrostatic potential and density perturbations. 
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In Sect. 3.2 we calculate the partial and total non-linear 
currents in terms of a given spectrum of irregularities. In 
Sect. 3.3, we calculate partial energy inputs and turbu- 
lent heating of electrons and ions and verify that general 
Eqs. (10) and (20) remain in this approximation exactly 
valid. 

3.1. First-Order Linear Wave Relations 

Fejer et al. [1984] studied the linear theory of collisional 
waves for strongly magnetized electrons and arbitrarily mag- 
netized ions. On the other hand, Buchert et al. [2006] as- 
sumed arbitrarily magnetized electrons but unmagnetized 
ions. Since no one has published the general 3-D linear rela- 
tions that cover all cases, we do this here and a more general 
version in the appendix of the companion paper [Dimant and 
Oppenheim, 2011]. In reasonable agreement with our PIC 
simulations, we assume that most of developed turbulence 
lies in the long-wavelength, low-frequency range of kU >C 1 
and oj <gC Vi, where uj is the wave frequency and U is the 
mean free path of ions with respect to dominant ion-neutral 
collisions. This allows us to employ a two-fluid, as opposed 
to kinetic, model and order various terms in the momentum 
equations with respect to the small parameters kh and u/vi. 
In this ordering, particle inertia and pressure gradients are 
second-order effects and can be neglected. This first-order 
approximation yields a dispersion relation and the corre- 
sponding relation between fluctuations of the plasma den- 
sity and electrostatic potential. These relations are common 
for all E / D-region plasma instabilities [Dimant and Oppen- 
heim, 2004]. The neglected second-order corrections are cru- 
cial for the linear wave dissipation and instability driving, 
but they are of less importance to the spatially/temporally 
averaged energy transfer and plasma heating. To the second- 
order accuracy, the general linear wave theory for arbitrarily 
magnetized plasmas is developed in the appendix of Dimant 
and Oppenheim [2011]. 

Though we discuss here the linear wave relationships, 
that does not mean that we consider the linear stage of 
instability. On the contrary, we assume a fully developed 
and non-linearly saturated turbulence in which the linear 
wave growth is balanced by non-linearities. However, we 
presume that these non-linearities only weakly modify the 
linear wave relationships, so that we include all non-linear 
terms, along with the instability driving or damping terms, 
into the second-order corrections and neglect their feedback 
on the first-order relations. 

Under the first-order approximation, each group of par- 
ticles has only two balancing forces: the Lorentz force and 
resistive collisional friction, 



q s (E + V 3 x B) = m s v s V s , 



(21) 



where s — e,i; we presume V s lies a neutral frame of refer- 
ence and an undisturbed magnetic field, B = Bo. Introduc- 
ing the conventional magnetization parameters, n s = Q s /i/ 3 , 
where Q s = \q s \B/m s are the electron and ion gyrofrequen- 
cies, we obtain from Eq. (21) 



V,, 



B 



B 



(22) 



K e [-E ± + K e (EX b)} - Ki[E ± +Ki(EX b)} 

V e± = 7^—, i Vil- = 



(1 + K*)B 



(l + /s?)B 



(23) 



with the subscripts ||,_L denoting the components parallel 
and perpendicular to B, respectively; here B = \B\ and 
b = B/B. 

Further, we calculate the relative mean velocity between 
electrons and ions, U = V e — Vi, which plays an important 



role in many relations. Equations (22) and (23) give 
U\\ = ^ , 



(24a) 



ft _ (K e + Kj) [(K e - Kj) (E X b) - (1 + KjKe) E±] . . 

(H-/s?)(1 + k?)B • (Mb) 
Reversing Eq. (24), we obtain 

Ei\ Uu 



B K, e -\- Ki 

E± (lie - Kj) (U X S) + (1 + KiK e ) U± 



B 



fci ~\~ Ke 



(25) 



so that Eqs. (22) and (23) yield the following mean drift 
velocities of electrons and ions in terms of their relative ve- 
locity: 



Ve, = 



K e [U± - Ki(U X 6)] 



Vii = - 



Kj[U± + K e (U X S)] 



(26) 

(27a) 
(27b) 



All these linear expressions can be applied separately to the 
undisturbed quantities and wave perturbations. 

Now we consider wave perturbations of the field, plasma 
density, and fluid velocities and apply to them the discrete 
Fourier transforms introduced by Eq. (11). Note that the 
neglect of non-linear terms with respect to wave perturba- 
tions leads to a unique fc-dependence of the wave frequency, 
u> = wg(fc), via the corresponding first-order linear disper- 
sion relation. Such unique dependence may break the pre- 
sumed discreteness of the wave frequency, u = 2irm/T. 
However, the quasi-period T can always be chosen suffi- 
ciently long so that the interval between adjacent discrete 
frequencies, Acu — 2n/T, becomes much less than the finite 
spectral width around u — cjg; see also the discussion below 
in the paragraph following Eq. (40). 

Expressing velocity perturbations, 5U ei g^, in terms of 
the electrostatic potential perturbations, via SE^ u = 

—ikSQ^^, we obtain from Eqs. (22) and (23): 



Ke[-k± + Ke(k± X &)]<?$£ ; 



(1 + K*)B 

Ki[k± + K t (k± x b)]6^ 
(1 -nf)B ' 



(28a) 
(28b) 

(28c) 



Writing the quasi-neutral continuity equations as 

d t n + V ■ (nVi) =0, V • (nU) = 0, (29) 

where we neglected wave variations of ionization-recombination 
balance, and linearizing them with respect to density per- 
turbations, Snj:^, we obtain 



no n ( J° 

k 



k-U 



(30) 
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Here Q.f = cj-fc-U* o is the linear wave frequency in the ion 
frame; 8U^ w = SV e g ^ - SV^ w , and, according to Eq. (24), 

(Ki + K e ) [(k e — Ki) (E X b) - (1 + KiK e ) E ] 



U ° (1 + k!) (1 + n\)B 

Using Eq. (28), we obtain 

k 2 



fc 2 



2^ S H» 



2 \ 



so that 



lK;K e (/t e + «») (1 + V>fc ) fe ± <^£, u 

(1 + K 2 )(1 + k 2 )B 



where 



v>± = 



I + (I + + K 2 ) 



V e Vi 



k\ 



(31) 

(32a) 
(32b) 

(33) 

(34a) 
(34b) 



The 2-D parameter ip± is conventional. The newly defined 
3-D parameter tpt generalizes the traditional 3-D parameter 
tp = ip±(l + «f fen / originally introduced for magnetized 
electrons, k 2 3> 1, and unmagnetized ions, s- « 1 [e.g., 
Farley, 1996]. In a more general case of k 2 ^> 1 but arbi- 
trary Ki, the parameter replaces the product (1 + K 2 )tp 
[Fejer et al., 1984]. In the lower ionosphere, the difference 
between the two is negligible, (1 + K 2 )ip — ipfe = K i V'x — 



6 = m e v e / (rtiiVi) 



x 10 [Dimant and Milikh, 2003; 



Dimant and Oppenheim, 2004]. 

Using Eqs. (32) and (33), we obtain from Eq. (30) the 
first-order relation between the linear fluctuations of the 
density and electrostatic potential, 



(1 + k 2 e )(1 + K 2 )B(k ■ U ) (8n % l 



KiK e (n e + Ki) (1 + 1p%)k^ 



no 



(35) 



as well as the expression for the first-order wave frequency 
in the ion frame, 

(i) = (l + K 2 )[l + (l + K 2 )fc 2 /fci](fc-i7o) 

E K e (« e + Ki) (1 + ^) ' ( ' 

In the neutral frame, the corresponding frequency, cjg = 
nf + k- V i0 , is 

(k b — m)(k • Uo) KiK e k ■ (Uo x b) . . 



(k c + Ki)(l + 1p%) 



K e -\- Ki 



When obtaining these expressions, we used the easily de- 
rived relation 

(1 + KiK e ) k± + (1 + kI)(1 + K 2 )k 2 = KiK e k 2 ± (l + l/j^). (38) 

As might be expected for arbitrarily magnetized electrons 
and ions, Eqs. (35) and (37) are symmetric with respect 
to the interchange between electrons and ions that requires 
Ki, e <r¥ K e ,i and Uo o — Uo- Equation (36) for the Doppler- 



shifted frequency in the ion frame, fig , is symmetric with 
respect to a similar expression for the wave frequency in the 



electron frame, fil e ' = cur — k ■ V e o = fi^ — k ■ Uo, 

k I* k 



Q (e) = _ (l + ^)[l + (l + « 2 )fc 2 /fc 2 ](fc-i7 C ,) 
k Ki (Ki + K e ) (1 + Vg) 



(39) 



Notice that for prevalent waves with positive k-Uo the shifted 
electron frequency fig 8 ' is always negative, while the corre- 
sponding ion frequency fil 1 ' is positive. This reflects the fact 
that such waves, regardless of the particle magnetization, al- 
ways lag behind the streaming electrons but move ahead of 
the ions. 

Despite the formal symmetry in the general relations be- 
tween electrons and ions, their actual contributions are not 
equivalent. In the lower ionosphere, since im ~ 30 amu 
and h> e /vi — 10 [Kelley, 2009], the ratio K e /Ki is huge, 
K e / Ki ~ 5500. Then practically everywhere, except for the 
D-region altitudes below 80 km, electrons are strongly mag- 
netized, K e ^> 1, whilst in most of the E / D-region electro- 
jet ions are largely unmagnetized, k, C 1, reaching only a 
partial magnetization, m > 1, above 110 km. Apart from 
the D-region altitudes but practically throughout the en- 
tire electrojet, setting K e S> 1 with arbitrary m, we reduce 
Eqs. (34) to (37) to simpler relations [Fejer et al, 1984]: 



1 + ri + K 2 )^h. 
1 + (1 + Kj ^fc 2 



k-Up 
1 + Vs 



— Kik ■ (Uo x b), 



(40) 



= i 



■mjVi(\ + K 2 j)(k ■ Up) 
efe 2 (l + ^) 



Sn s 



k,u> 



no 



For unmagnetized ions, Ki <C 1, Eq. (40) reduces further to 
the conventional relations [Farley, 1996]. 

Recall that all these relations represent only the first- 
order approximation with respect to the small parameters 
kU and u>/vi. Neglect of pressure gradients, particle inertia, 
etc., results in the highest-order dispersion relation for the 
real part of the linear wave frequency, Eq. (37). It contains 
no 7r/2-phase shifted corrections that determine the linear 
wave growth or damping. For arbitrary particle magnetiza- 
tion, the corresponding second-order corrections to the dis- 
persion relation for the FB and gradient drift instabilities are 
obtained in the appendix of Dimant and Oppenheim [2011]. 
We should bear in mind, however, that our quasi-periodic 
description of the non-linearly saturated steady state re- 
quires all wave frequencies to be strictly real. This implies 
that the second-order destabilizing linear factors are on av- 
erage balanced by non-linearities. Furthermore, non-linearly 
saturated turbulence itself has a strongly time- varying char- 
acter caused by dynamic mode-coupling [e.g., Hamza and 
St.-Maurice, 1993a, b; Dimant, 2000]. All this breaks the 
narrowband approximation for which any wave frequency u, 
via the linear dispersion relation, is uniquely determined by 
the corresponding wavevector k, uj w uj%. This means that 
the (5-function dependence of the frequency spectrum is ac- 
tually spread over a finite band around cjg. Accurate theo- 
retical description of such non-linearly saturated states is ex- 
tremely difficult, especially for strong turbulence generated 
by the driving field well above the instability threshold. To 
avoid serious difficulties and mathematical complexities, we 
will continue using first-order linear relations like Eqs. (28) 
and (35) and neglect next-order linear and non-linear correc- 
tions. We expect the corresponding errors to be reasonably 
small, although these expectations need additional testing. 

3.2. Non-linear Currents 

As shown in Sect. 2, the non-linear current (NC) plays 
an important role in energy conversion. To further clarify 
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its physical meaning, we calculate separately NCs for elec- 
trons and ions, J^ c (s = e,i), in terms of a given density- 
irregularity spectrum, 8n^ w . These partial currents will be 
used in Sect. 3.3 for comparison of the partial turbulent 



Joule heating, Eo ■ , with the corresponding frictional 
heating terms. 

Using Eqs. (12), (28), and (35), we obtain the partial elec- 



tron NC density, 

(l + «e) 



Je = 



-e{8nSV e ) 



1 + Ki ) en 



Ki (^-e ~\~ Ki) 



E 



^0 



k-U 



fell + 



k± — K e (k± x b) 



1 + k1 



no 



(41) 



7NC 



Similarly, we obtain the partial ion NC density, j t 
e{5nSVi) = e^^^fej^j, which differs from Eq. (41) 
by the replacement K e ,i <r¥ Ki, e and the "plus" sign in 
front of (jfcx x i). The total NC, j^ = -e(5n5U) = 

- e J2k,^o Sn *k SU k' is then s iven b y 



J = Je + Ji = 



en 



5n k,u 



(1 + k\) (1 + Ki) fc|| + (1 + K,K e ) fcj 



n 

- {He 



■ Ki) (k x b) 



(i + ^)fei 



(42) 



At altitudes above 100 km, where both n e ^> 1 > Ki and 
KiK e — ipj 1 >• 1 hold together, according to Eq. (31), we 
have 

Box!) - KiEo 



U = 



(43) 



while the ion non-linear current turns out to be negligible 
compared to that of electrons, J^ c « J^ . We expect no 
spectral asymmetry along B, so that w - /(fcj_)fc|| = 0. 
As a result, the NC density, J* 10 _L B, reduces in this limit 
to 



en 



E 

k,u;^0 



(fc x b)(k-U ) 



n 



(44) 



Within the bulk electrojet where ions are unmagnetized, 
k, < 1, the fc-spectrum of irregularities is largely perpendic- 
ular to Eo, so that j 1 * has there a predominantly Pedersen 
direction. This is of paramount importance for the global 
MI coupling [Dimant and Oppenheim, 2011]. 

3.3. Partial Energy Inputs and Turbulent Heating 

Calculating the partial energy inputs and frictional heat- 
ing sources for a specific plasma model and given turbulence 
spectrum allows to quantitatively understand how the tur- 
bulent energy is distributed between different plasma com- 
ponents. Here we obtain such expressions for arbitrarily 
magnetized two-fluid plasmas in the quasilinear approxima- 
tion and verify, in particular, that exact Eqs. (10) and 
(20) remain exactly valid. The turbulent frictional heat- 
ing sources found here could be included into ionosphere- 
thermosphere computer models, as explained in the com- 
panion paper [Dimant and Oppenheim, 2011]. 

Multiplying Eq. (21) by the plasma density and corre- 
sponding fluid velocities gives 



E ■ j s = m s u s nV s 



(45) 



(s = e, i), where j e = —enV e and ji = enVi are the electron 
and ion partial current densities. These expressions relate 
the work done by the field E on the s-particle currents to the 
corresponding sources of frictional heating. As mentioned in 
Sect. 2, the actual frictional heating of s-type particles is 
smaller than m s v sn nV 2 because a fraction of the acquired 
field energy equal to m s /(m s + ra„) goes immediately to 
the colliding neutrals without heating the plasma particles 
[Schunk and Nagy, 2009]. This is especially true for ions 
with mi ~ m n , whose frictional heating source is nearly half 
of miUinVi . In what follows, we will refer to m a v an nV 2 
as the s-n 'heating term' with the caveat that the actual s- 
particle frictional heating is described by m n /(m s + m„) of 
m s u a nV s 2 , while the remaining fraction goes to neutral (n) 
frictional heating [Dimant and Oppenheim, 2011]. 

Equation (45) includes total field energy losses and 
plasma heating. The zero-order heating alone, i.e., that 
without any plasma turbulence, is described by 



E ■ jso = m s v s noV sQ . 



(46) 



Subtracting Eq. (46) from Eq. (45) and averaging yields 
the turbulent heating rates, 

Eo ■ jT C + (5E ■ Sj s ) 

= m a v s (n {SV 2 ) + 2V s0 ■ (SVJn) +{SV s 2 Sn)), (47) 

where J^ c = —e{5nSV e ) and jf c = e{8nSVi) are the par- 
tial electron and ion contributions to the total NC given by 
Eq. (42). The LHS of Eq. (47) describes the additional 
average work by electric fields in the turbulent plasma on a 
given plasma species. The RHS represents the correspond- 
ing average turbulent heating per unit volume. It is not 
obvious, however, that the approximate quasilinear expres- 
sions obtained above ensure that the two sides of Eq. (47) 
are really equal. To verify the equality and clarify the physi- 
cal meaning of various terms, we will calculate the two sides 
of Eq. (47) separately. 

We start by calculating the RHS of Eq. (47) that de- 
scribes the frictional heating. The positively determined 
and dominant term m a v s no{8V 2 ) describes energization of 
individual plasma particles, while the other terms are asso- 
ciated with density variations at a given location. Consis- 
tency of the quasilinear approximation requires neglecting 
the last, cubically non-linear, term. To calculate average 
quadratically non-linear quantities in terms of given irregu- 
larity spectra, we use Eq. (12). 

By using Eq. (28), for the dominant electron and ion 
heating terms, we obtain 



m s u a no (SV 2 ) = 



enoK a 



E 



fc 



1 + K§ 



(48) 

If k e 3> Ki and K e ^> 1 then the ion perpendicular heat- 
ing dominates over the corresponding electron one, while 
for the parallel heating the reverse is true. At lower D- 
region altitudes, where electrons are partially demagnetized, 
Ki <JC K e < 1, the electron heating always prevails. Using 
Eq. (35), we rewrite Eq. (48) in terms of a given density- 
irregularity spectrum as 



e v e no (SV 2 ) 



en B(l + K 2 e )(l + k 2 ) 2 

(/v e -f- Ki^K^Ke 

<l + (l + K 2 e )k 2 /k 2 ± ](k-U ) 2 



x £ 



no I 
(49) 



and similar for ions with the symmetric replacement K e ,j o 
Ki te . Using Eqs. (35) and (38), for the combined heating 
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rate we obtain 

m e u e n (<5V e 2 ) + m,wio (SV?) 

eBn ip±{l + k 2 )(1 + k 2 ) (fc ■ Up) 2 



no 



(50a) 



ena (m + k c 1 - - - 



rSm) Ea+^M - (sob) 



BVx(l + ^)(l + ^) 

Equation (50), presented here in terms of both Sn% ^ 
and <53>£ , allows one to understand the quantitative dif- 
ference between 2-D and 3-D turbulent heating. The 3- 
D effect of parallel-field dominated turbulent heating [Di- 
mant and Milikh, 2003; Milikh and Dimant, 2003] is en- 
tirely due to the difference between tp^ and ip± in the mul- 
tipliers (I + tpt)- While the total perpendicular turbulent 
heating is described by the (I + ip±) component, the to- 
tal parallel turbulent heating derives from the remainder, 
ip^-ip^ =i/>x(l + Ke)(l + «i )fcj| /k± , as seen from Eq. (34) . 
If |5ng^| 2 were equal in 2-D and 3-D then, according to 
Eq. (50a), larger (1 + ipA would reduce total turbulent 
heating in 3-D compared to 2-D. We have, however, the re- 
verse by the following reason. Heuristic arguments require 
the perpendicular rms turbulent fields ~ U , Q k\\8§^ J 2 
to be approximately equal in 2-D and 3-D [Dimant and 
Milikh, 2003; Dimant and Oppenheim, 2011]. According 
to Eq. (50b), this leads to stronger heating in 3-D com- 
pared to 2-D, entirely due to larger density perturbations, 
E fe W„ \Sn s J 2 oc £ Wo (l + ^mi^.J 2 - The same 
pertains to the non-linear currents discussed in Sect. 3.2. 
All these heuristic inferences have been confirmed by our 
recent supercomputer PIC simulations [Oppenheim et al, 
2011]. 

The second term in the RHS of Eq. (47), calculated 
by expressing 8V e , i and V e ,io in terms of SU and Uo, us- 
ing Eq. (27), and Uo x b in terms of E , using Eq. (25), 
is 



2m,i ( ViVio ■ {SnSVi) ) 
= 2en (l + K 2 e ) 

f^e ( ~4~ f^e ) 



E 

k,ui^0 



(k-E )(k-U ) 



5n k,u, 



n 



(51) 



and similar for electrons, with the symmetric replacement 
At higher altitudes, where tp± — (n e m) 1 < 1 
and k e > 1, the ion component dominates over the electron 
component, while at lower altitudes, ip± > 1, the reverse 
holds. Adding up the two plasma components gives 



2m i v i (y i o ■ {SnSVi)) + 2m e v e (V e o ■ {8n8V e )) 



■ 2eno (1 + ipi 



■ E ■ 



k-Eo)(k-U ) 



Sn k,» 



no 



(52) 



In the most of electrojet, these terms are small compared to 
those in Eq. (50), although at the top of the electrojet they 
can become comparable. 

To verify that general Eqs. (10) and (20) exactly hold in 
our quasilinear calculations, now we proceed to calculating 
the LHS of Eq. (47). Each of the electron and ion terms, 
{SE ■ 5j e ,i), can be separated into two distinct parts, 



{SE ■ Sj a ) = {5E ■ 6j.)i + (SE ■ Sj.) 2 , 



(53) 



where {SE ■ 6j 3 )i = -en {SV s ■ SE), {SE ■ Sj 3 } 2 = -eV e0 • 
{SnSE) , and we neglect the cubically non-linear corrections, 
e{SnSV a -8E). Using Eqs. (25) and (28), we find that the first 



term, {SE ■ 5j s )i, equals the fixed-density frictional heating 
rate, 

{6M-63.) 1 = m t v.n (5V?), (54) 

while the combination of the electron and ion second terms 
in Eq. (53) yields 



{SE ■ 8j e ) 2 + {SE ■ 5ji)2 = {SE ■ 8j) 2 
eBn (1 + k 2 ) (1 + k\) ip± 



E 

k,ui^0 



(fc ■ Uof 



n 



(55) 



By combining Eqs. (50) and (53)-(55), we verify that, in 
accord with Eq. (10), {SE ■ 8j) — 0. Due to this, turbulent 
Joule heating, i turb = {E-j)- E ■]o = E -] NC + {SE ■ Sj), 
in accord with Eq. (20), equals the total average work of 
the external electric field on the total non-linear current, 
P N c = Eo ' J^ - This quantity, in turn, should equal the 
total frictional heating source given by combining the RHS 
of Eq. (47) for all plasma species. Indeed, using Eq. (42) 
for j^ , we obtain 



Pnc 



eBn ip±(l + k 2 )(1 + k 2 
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+ 2(1 + V±)en 
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k,ui^0 



(l+^)fcl | no 



(56) 



which equals the total turbulent heating found by combining 
Eqs. (50) and (52). 

Thus, our fluid-model calculations for arbitrarily magne- 
tized plasma particles in the quasilinear approximation have 
fully confirmed the exact general relations derived in Sect. 2 
exclusively from Maxwell's equations. The partial relations 
for various energy conversion terms show how the deposited 
field energy is divided between the electron and ion heating 
channels and allow one to properly interpret the correspond- 
ing terms. 

4. Global Energy Flow 

Now we discuss how the energy from the Earth's mag- 
netosphere is deposited to the i?/D-region electrojet. From 
the above treatment it is clear that the global energy flow 
should be analyzed based on the total average currents that 
include the NCs. To obtain energy balance equations similar 
to Eqs. (2) for spatially and temporally averaged fields in- 
troduced in Eq. (7), we use separate linear Maxwell's equa- 
tions for the average fields and for disturbances caused by 
the total average currents, (j) = j tot = jo+jNL, as we did in 
deriving Eq. (9) for the wave perturbations. As a result, ig- 
noring the induction component of (E) = Eq, i-e., assuming 
V r - x Eo — 0, we obtain 



£0 (£ 2 + cW 



+ V r - • AS = -Eo ■ j 



(57) 



where r, rare the large-scale variables introduced in Sect. 2, 
and 

AS = Eo X AB = e c 2 E x AB. (58) 
Mo 

Here AB = {B) — Bo is the quasi-stationary large-scale mag- 
netic field disturbance, where by Bo we mean the geomag- 
netic field undisturbed by the electrojets, V? x Bo = 0, so 
that Vp- (Eo x B ^j = 0. The disturbance AB is caused by 
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the total average electrojet currents, V x AB « noj tot - It is 
usually so small, |AB| < 10~ 3 Bo, that can be neglected in 
the expression for the magnetic energy density, (B) 2 « Bq. 
However, AB is crucial for the Poynting flux AS* that pro- 
vides a downward flow of field energy from the magneto- 
sphere to the E j £ )-reg ion ionosphere. Further, bearing in 

mind that j tot = <r tot • Eo, the work of the external field 
on the total current in the RHS of Eq. (57) can be writ- 
ten as E ■ j tot = <Tp ot £o, where and a P ot are the total 
conductivity tensor and Pedersen conductivity, respectively 
[Dimant and Oppenheim, 2011]. In the turbulent electrojet, 
the quantity ap ot E 2 represents the combined laminar and 
turbulent Joule heating. Bearing in mind quasi-stationary 
conditions, from Eq. (57) we estimate the magnitude of the 
Poynting flux on top of the electrojet as \AS\ ~ a P ot _EoI/||, 
where Ly is the characteristic size of the electrojet along the 
nearly vertical magnetic field. The corresponding distur- 
bance of the magnetic field, AB, has a significant compo- 
nent perpendicular to both Eq and Bo with the magnitude 

~ /ioOp ot i5o-I'|| • 

Now we do the same for the plasma. Averaging Eq. (5) 
over small turbulent scales and neglecting (SE-8j), in accord 
with the analysis of Sect. 2, we obtain 

d T + v - ■ Z)<*.> = £° ■ f ot + E< L *>- ( 59 ) 

s s s 

While the RHS of Eq. (57) has only one term describing 
the total field energy loss to charged particles, the RHS of 
Eq. (59) has two terms. The first term, E ■ j tot = a^E 2 , 
describes the total Joule heating rate due to the total energy 
input from the fields. The second term describes particle en- 
ergy dissipation in the abundant neutral atmosphere. For 
particle fluxes dominated by the thermal bulk, the two terms 
almost cancel each other locally, so that the flux divergence 
term V r • ~}2 S {K S ) is expected to be much less than Vf AS 
in Eq. (57). This means that the energy transport of iono- 
spheric particles should be much weaker than that of fields. 
The reason why the particle energy transport plays a minor 
role is associated with relatively short mean electron and 
ion free paths, even in the almost vertical direction along 
Bo- We should bear in mind, however, that the fluxes (K s ) 
can also include high-energy precipitating particles that may 
comprise a significant part of the field-aligned (Birkeland) 
currents within auroral arcs. These occasions may break the 
nearly perfect local balance between the Joule heating and 
atmospheric cooling. Although in these cases |V r • ^2 S (K S )\ 
should still remain less than AS\, the two fluxes may be 
closer to each other than those for the nearly local thermal- 
bulk dominated balance. 

Now we compare the total field and particle fluxes along 
the nearly vertical magnetic field in greater detail. We can 
estimate the particle energy flux magnitude as | ^ s K s \ ~ 
{£)j\\ , where {£ ) is an effective energy of charged parti- 
cles. Using the charge flow conservation, V • j = 0, we 
estimate the parallel current density as jy ~ L^apEo/Lp, 
where Lp is the characteristic scale of current density vari- 
ations in the Pedersen direction. As a result, we have 
IE 3 (^)I/I A <S| ~ (£)/(e|A$|), where A$ ~ L P E is the 
characteristic cross-polar cap potential. Typical values of 
A$ are ~ 100 kV. If this current is mainly provided by 
thermal bulk particles then their energy is many orders of 
magnitude smaller, {£) < 1 eV. Precipitating particles can 
have a rather high energy, (£) < 30 keV [e.g., Ashrafi et al, 
2005], so that the particle energy flux can sometimes be 
comparable to the Poynting flux. 

Bucheri et al. [2006] in their analysis of global energy flow 
between the magnetosphere and £-region ionosphere put 



more stress on the Birkeland currents. These field-aligned 
currents are important for charge conservation and they pro- 
vide the MI coupling via particle precipitation. Further- 
more, electron fluxes along Bo, because of their high parallel 
mobility, ensure effective mapping of the electric field from 
the magnetosphere to lower ionosphere. From the energy 
transport viewpoint, however, our analysis gives preference 
to the Poynting flux. It is interesting to note that 2-D or 
3-D simulations of the E / D-region instabilities in purely pe- 
riodic boxes have always generated the DC currents j tot in 
a plane perpendicular to Bo- According to Maxwell's equa- 
tions, j tot must generate quasi-stationary loop-like magnetic 
disturbances AB, which in the employed electrostatic codes 
are ignored by definition. Furthermore, such AB would even 
violate the imposed periodicity, unless one presumes on the 
3-D box boundaries surface currents that exactly balance 
the volumetric currents. If, however, one did not ignore AB 
then the corresponding Poynting flux, c 2 eoEo x AB, would 
be directed inward the box and provide the required energy 
input. 

5. Summary and Conclusions 

Plasma turbulence generated by E/D -region instabilities 
can contribute significantly to the global energy exchange 
between the magnetosphere and ionosphere. The spatially- 
temporally averaged energy deposit in the turbulent electro- 
jet is given by Eo ■ j tot , where the total current density j tot , 
in addition to the regular current density, jo = ffp • Eo, in- 
cludes the non-linear current (NC) density, J^ , caused by 
low-altitude plasma turbulence [Rogister and Jamin, 1975; 
Oppenheim, 1997]. The work of the external field on the 
NC, Eo • j^ , provides the required energy input for anoma- 
lous (turbulent) heating of both electrons and ions [Bucheri 
et al., 2006]. In Sect. 2 we prove this earlier conjecture using 
first principles with virtually no approximations. Specific 
fluid-model calculations for arbitrarily magnetized plasma 
in the quasilinear approximation show how exactly the de- 
posited field energy is distributed between electrons and 
ions. This yields explicit NC expressions given by Eq. (42), 
and turbulent sources of anomalous electron and ion heating 
given by Eq. (49). It has been known for a long time that 
anomalous electron heating (AEH) is largely a 3-D effect 
caused by turbulent electric fields parallel to B [St. -Maurice 
and Laher, 1985; Providakes et al., 1988; Dimant and Mi- 
likh, 2003; Bahcivan et al., 2006]. According to the new 
development, this requires significantly larger Eo ■ J 1 " 10 in 
3-D than in 2-D. Due to the mirror symmetry along B, nei- 
ther Eo nor ~j NC have any significant components in the B 
direction. The difference in Eo ■ J^ between the 3-D and 
2-D is entirely due to noticeably larger 3-D NC caused by 
larger density perturbations, as explained in the text below 
Eq. (50). This prediction has been confirmed by our recent 
2-D and 3-D PIC simulations [Oppenheim et al, 2011]. As 
discussed in the companion paper [Dimant and Oppenheim, 
2011], a strong NC, either by itself or combined with AEH 
[Milikh and Dimant, 2003; Milikh et al, 2006], can signif- 
icantly increase the global ionospheric conductances. This 
may have serious implications for predictive modeling of MI 
coupling and space weather. 

Appendix: Validation of Electrostatic 
Appoximation 

Using some relations obtained in Sect. 2, we now discuss 
the validity of the electrostatic approximation in describing 
wave processes in the lower ionosphere. To the best of our 
knowledge, the electrostatic approximation for the B-region 
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instabilities has always been employed but no detailed anal- 
ysis of its validity been done, especially for the general span 
of magnetization conditions considered in this paper. There 
are occasions when induction fields in the ionosphere play 
an important role [e.g., Amm et al, 2008; Vanhamaki et al., 
2007]. The main reason why one might question the va- 
lidity of electrostatic approximation in the description of 
S/-D-region wave processes is that even tiny induction cor- 
rections to the wave electric field along B could modify the 
small parallel turbulent electric field largely responsible for 
AEH. 

To estimate a possible non-clcctrostatic component of a 
separate wave field harmonic, SE^ , we split it into the elec- 
trostatic (curl-free) part, 8E^ S ^ = (k-SE^ w )fc/fe 2 , and the in- 



duction part, 5E l fV = 
Then Eq. (15) yields 



k.u 



k.u 



SEf = 

k,u 



SEf D = 

k.U 



»( fc • 5 k,J k 

£0Ujk 2 

iiok x (k x Sj^ t 



£o(fe 2 



-kxikxSE^J/k 2 . 

(60a) 
(60b) 



)k 2 



8E t 



The wave current density, 8jt ^ , is determined by an 
anisotropic response to the total turbulent electric field, 
SE^ . Even if the turbulent field is largely electrostatic, 
SE^ || k, this anisotropy can give rise to a non- 
negligible component of Sjj; w which is not parallel to k. Ac- 
cording to Eq. (60b), this component generates non-zero 
k x (fe x 5j %w ) = k(k ■ J - k?6j^ u , i.e., finite SE^°. 

Using a perturbation technique in which the zero-order 
wave field is electrostatic, 5Er,, « SEr S = —ikS^s,,, 
one can easily estimate the next-order induction compo- 
nent, 8E™^ . Considering the wave with k in the x, 2-plane 

(z||Bo), neglecting density perturbations, and applying to 
5 E% u the regular linear conductivity, we obtain 



where we have neglected op compared to erf? and k« com- 
pared to k\ w k 2 . We also have taken into account that 
in E/D-region processes the wave phase speed, V p h = u/k, 
is many orders of magnitude less than the speed of light, 
so that d 2 in the denominator of Eq. (60b) can always be 
neglected compared to k 2 c 2 . 

Comparing Eqs. (62) and (63) for both parallel and per- 
pendicular to Bo components and presuming crjjfcy/fex < cr P , 
wo SCO that the smallness of \6E5* D /5Eg S | is ensured by 



Vr 



P h 



1/2 



Vr 



P h 



1/2 



(64) 



The first inequality pertains to the parallel components of 
the turbulent electric field, SEr .. , while the second, less 
restrictive, inequality pertains to the corresponding perpen- 
dicular components, SE^ ± . As mentioned above, break- 
ing of the parallel-field condition may be of importance be- 
cause 8E^,, is crucial for AEH. For magnetized electrons, 
up to the top electrojet where K e 3> 1 and m ~ 1, we have 
tTp ~ Kine/B, (j^ ~ ne/B, crjj « K e ne/B, so that the most 
restrictive first inequality becomes 



^«f-) 1/2 ^e 

c UJ 



1.4 x 10" 



(65) 



Even for an extremely strong convection electric field of 
Eq ~ 150 mV/m, the wave phase velocity may reach 3 km/s 
at most, Vph/c ~ 10 -5 , so that this condition holds with a 
big reserve. Thus for typical E/D-region wave processes the 
electrostatic approximation is always valid even for small 
parallel electric-field components, SEg ... 
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Since the anomalous conductivity discussed in the compan- 
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